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A method is presented to generalize the power detectors for short bursts of gravitational waves 
that have been developed for single interferometers so that they can optimally process data from a 
network of interferometers. The performances of this method for the estimation of the position of 
the source are studied using numerical simulations. 
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. I. INTRODUCTION 

5-H ; 

^-i. The most interesting advances in gravitational wave (GW) astronomy are likely to come from the cooperative 

' operation of the interferometric detectors of gravitational radiation that are currently in an advanced stage of com- 

. missioning in North America, Europe, and Asia. The two principal advantages of operating the interferometers in 

' a network are the significant increase in observational sensitivity, and the availability of position information for the 
source of the gravitational radiation. As discussed in a realistic goal for Advanced detectors is to trigger the 

' observation of the electromagnetic radiation associated with compact binary mergers. 

] The simplest approach to analyzing data from a network of GW detectors is to process the data from each inter- 

^ I ■ ferometer individually, and to combine the results incoherently. Event lists are generated for every interferometer, 

I I and coincidence gates of various levels of complexity (which may include information on time, frequency, amplitude, 

j- . etc.) are used to combine the lists into a single set of GW candidates. The source position can be estimated by 

' triangulation when good information about the arrival time of the bursts in all interferometers is available. This 

CO incoherent approach has been used for all the analyzes performed on real data to date 0, 1^ • 

Recent theoretical work has demonstrated the existence of an optimal statistic for a coherent analysis when matched 

Q \ filtering is possible 0,0, i.e., when the phase evolution of the source is precisely known. In essence, the coherent 

O^" analysis combines the data before the search for the presence of a gravitational wave is performed. This paper describes 

I I a parallel approach for the more complex case where no precise model of the gravitational wave signal is available. 

bJ[)' Such signals are often detected in single interferometers using power detectors which rely on a quadratic form of the 

■ data as a detection statistics [1, 0, |M IS ^ • I describe in section ^1 an algorithm which optimally combines 

• i-H , data from a network of interferometers into a single, synthetic data stream which can be processed by any of these 

' single interferometer power detectors. Section Hill presents results for some special cases studied analytically or with 

^ \ numerical simulations. A more complete discussion of the algorithm and of its performances is presented elsewhere 

■ - - ' Il3, and is summarized in section Hvl 



II. ALGORITHM 



Consider one of the many power detectors that were developed to detect GW bursts of arbitrary waveform in the 
data from individual interferometers. It is assumed that the power detector produces a estimate P{y) of the power 
present in the signal for a given data vector y, and that this estimated power can be expressed in terms of a scalar 
product, 

P{v) = \v? = vy- (1) 
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One may try to recycle such a power detector for a network of N independent interferometers by feeding it with a 
synthetic data vector Y , formed by hnearly combining the time-shifted data from the N detectors: 

N 

y = ^a,T(<5,)y,. (2) 
1=1 

Here, T is the time-shift operator, a is a set of scalar weights, and (5 is a set of time-shifts. 

The cross-correlation function between the plus and the cross polarizations of the GW signal (s_|_ and Sx) is defined 

by 

i?+x(A) =s+.T(A)sx. (3) 

Under the condition that this function has an extremum at A = 0, i.e. that its derivative satisfies (A)|a=o = 0, 
it can be shown that the signal-to-noise ratio (SNR) for the network power, P{Y), is maximum when the time-shifts 
S are chosen to correspond to the true position of the source (as expected), and when the weights a are a solution of 
the following eigenvalue problem 

N 

Xa^af + aAF+F+s+ ■ s+ + {F+F^' + Fl'F+)s+ ■ Sy, + F.^'^F^'^Sy, • Sx] = 0, (4) 

where i — 1,...,N. Here, A is an eigenvalue, af is the variance of yi, and F^ and F^^ are the interferometer 
beam-pattern functions for the two polarizations for the i*'' detector. 

Consequently, when i?^^(A)|A=o — Oi the SNR can be maximized with respect to only four parameters: the 
position of the source (6', </>), the ratio of the RMS of the cross polarization to the RMS of the plus polarization 
(A^/x = l*+l/l*x I) and the amount of "overlap" between the two polarizations (A+.x = s_|_ ■ Sx/|s_|_||sx |). It is very 
significant that the maximization can be performed over only four parameters; without this collapse of the parameter 
space from a space of extremely large dimensionality where the waveforms for the two polarizations are fully specified 
to a space with only four dimensions, no practical implementation of the method would be possible. The coherent 
power filter (CPF) algorithm described in this paper can search over all four parameters; it estimates the source 
parameters by returning the point in parameter space which provides the largest value of the SNR, as estimated from 
the data. Explicitly, the CPF algorithm is defined as follow: 

1. Pick a set of trial source parameters 9, (p, A^/^, and A+.x- 

2. Compute it from Eq.Q), and set 6=0. 

3. Form the synthetic response Y using Eq.(|2l. 

4. Use a single interferometer power detector to calculate P = \Y\'^. 

5. Retain the source parameters 9, 0, A+/x: and A+.x if P is the largest to date. 

6. Go back to Step 1 unless all points in parameter space have already been visited. 

The fact that the problem can be reduced to a maximization over only four parameters can be understood by 
a simple geometrical picture. Representing the time-series as vectors in a space of high dimensionality, the data 
observed by every interferometer are the sum of a random vector and of the linear combination of two signal vectors 
(for the plus and cross polarizations), appropriately rotated to account for the time of arrival of the signal at te 
detector. The weights in this linear combination are just the beam-pattern functions for each interferometer. If one 
picks a set of time-shifts that cancels out perfectly the differences in time of arrival, all the signal components of the 
time-shifted observation vectors will lie in a single hyperplane spanned by the two signal vectors. The orientation of 
that hyperplane is unknown, since the signal vectors are unknown. However, the goal is to maximize the norm of the 
weighted sum of the signal components of the time-shifted observation vectors. This maximization can be performed, 
because this norm depends only on the beam-pattern functions, the angle between the two signal vectors (i.e., A+.x), 
and the ratio of their norms (i.e., A+/x)i up to an unimportant overall scale factor. 
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III. EXAMPLES 

The examples presented in this section are all for the network formed by the LIGO interferometer near Hanford, 
Washington, the LIGO interferometer near Livingston, Louisiana, and the Virgo interferometer (the HLV network). 
Results for the study of a long, nearly monochromatic signal similar to a chirp from a binary inspiral are first 
discussed in relation to possible systematic errors on the position estimations performed by CPF. Results from a 
number of numerical simulations with a short, broader-band signal are then presented, in order to get insights into 
the statistics of the CPF algorithm as a position estimator. 

Obviously, the condition (^)|a=o = 0, which is assumed by the CPF algorithm, cannot be expected to hold true 
in general, for an arbitrary GW signal. I consider for instance a signal with characteristics similar to those of the chirp 
from a binary inspiral. It is a long, monochromatic signal with a 40 Hz frequency and a 7r/2 phase offset between the 
plus and the cross polarization waveforms, so that i?^x(^)|A=o ^ 0. A linearization of the normal equations for the 
SNR maximization problem can be used to show that the systematic error in position introduced by using the CPF 
algorithm for this signal is rather large, i.e. is above 10~^ radian for 99% of the sky. On the other hand, the size of 
the second order systematic error is negligible (< 10~^ radian) for approximately 25% of the sky. Since these errors 
are deterministic, a knowledge of the value of R'^^{A)\^^q is sufficient for implementing a correction procedure that 
removes the first order errors. This correction is feasible for ~ 25% of the sky, thus limiting the sky coverage of the 
CPF algorithm for source position estimation to one quarter of the sky, at any given time. For signals less sensitive 
to systematic errors than the one described in the example above, the noise in the detectors will obviously dominate 
the errors in the position estimates. Such errors are best studied using numerical simulations. 

In all numerical simulations, the noise in the three interferometers was taken to be white noise of unit variance, and 
the data were processed in 10 s long segments sampled at 16 384 Hz. The plus and cross polarization waveforms were 
independent realization of a bandlimited white noise process, of duration 1/16 s, and with negligible power below 125 
Hz and above 150 Hz. The plus polarization waveform was produced from a broadband noise of variance 2, and the 
cross polarization from a broadband noise of variance 1/2, so that on average A+/x = 2. Because of the independence 
of the two polarizations, on average A+.x was equal to zero. In all cases, the signal was rather strong; the signal was 
scaled so that its optimal filtering amplitude network signal-to-noise ratio^ was 13.4. All signals were injected from 
the northern hemisphere normal to the plane containing the three interferometers of the network. 

All simulations were performed using the LIGO Data Analysis System (LDAS) implementation of the CPF algorithm, 
which relies on the use of the TFCLUSTERS algorithm as the single interferometer power detector which measures the 
network power. In order to keep the size of the simulations at a manageable level, it was assumed that the value of 
A+/x was known exactly, and this parameter was not included in the maximization of the network power. For the 
other three parameters, the network power was first computed over 10 different values of A+.x uniformly distributed 
over the interval [—1, 1], 100 different values of the right ascension uniformly distributed over [0, 2tt[, and 100 different 
values of the sine of the declination uniformly distributed in [—1, 1]. Consequently, one parameter scan of the CPF 
algorithm involved 10^ runs of the TFCLUSTERS algorithm, and took ^ 225 s on a cluster of 32 Pentium IV computers 
clocked at 2 GHz and equipped with 512 Mb of RAM. In order to obtain a position estimate which was not limited 
by the resolution of the parameter space, a second run of the CPF algorithm was performed. For that second run, 
the value of A+.x was fixed to the value obtained in the first run, and both the right ascension and the sine of the 
declination were scanned with 50 points in each direction, in a square window of size 0.2 rad in both directions, 
centered on the position estimated in the first run. Figures ^ and |2 show an example of the two runs of CPF for a 
single simulated data segment. 

Running repeatedly the simulations gave a number of position estimates, which are plotted in figure|31 A significant 
number of points cluster around the position where the signal was injected, or around its mirror image with respect 
to the HLV plane, a natural consequence of the ambivalence present in any position estimation with only three 
interferometers. Approximately 25% of the trials had errors smaller than one degree, where the error is the arclength 
along the great circle connecting the estimate and the signal position or its mirror image, whichever is smaller. Most of 
the other trials gave positions that were significantly off the source position or its mirror image, but which lied on the 
equal-delay curve for the Hanford-Livingston interferometer pair. Physically, this results from the very good alignment 
of these two detectors, and the relatively poor alignment of the Virgo interferometer with them. In particular, it can 
be seen that in addition to hugging the Hanford-Livingston equal delay curve, the position estimates gather into 
clusters along that line. These clusters correspond to positions where the signals in the interferometers at Hanford 
and Livingston are in phase, and are at a phase offset with respect to the signal at Virgo which is an integer multiple 



i.e., the signal-to- noise ratio that could have been obtained has the arrival time, duration, and frequency band of the signal all been 
known beforehand 
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FIG. 1: An example of the output of the first pass of the CPF algorithm. The position plane with the value of A+.x giving 
the maximal network power is plotted. The gray scale represents the network power, normalized to its maximum value. The 
continuous curves show the loci of equal time delays for the three independent interferometer pairs. 
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FIG. 2: Same as Figure for the second pass of the CPF algorithm. The black dot shows the position with maximum network 
power. 

of the characteristic wavelength of the signal. 

IV. DISCUSSION 

There are significant difficulties inherent to the analysis and detection of GW bursts in the data from a realistic 
network of interferometers, such as the HLV network. The complications come essentially from the entanglement of 
the two unknown polarizations of the GW signal, and are exacerbated when the two polarization waveforms do not 
show significant differences in their amplitudes or structures, or when the interferometers are poorly aligned. A simple 
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FIG. 3: Scatter plot of the estimated source position, for 240 realizations of the simulation. 

and optimally efficient way of mitigating these problems was described in this paper. It can be used to efficiently 
recycle burst detectors that have been shown to be efficient at picking signals in a given data stream, and that are 
robust and well characterized when operated on non Gaussian or colored noise. 

Short bursts 10 cycles) have been studied numerically. For reasonable burst amplitudes, it was shown that their 
source could be located with good accuracy 1°) a significant fraction of the time. Variation of the parameters of 
the simulations |l2j | revealed that the position estimation error was scaling weakly with the burst amplitude, but was 
a stronger function of the degree of polarization of the wave, which was approximately controlled by the parameter 
A+/X. 

The misalignment of the Virgo interferometer with respect to the two LIGO detectors in the HLV network was also 
a significant source of error in the estimation of the position. In the simulations, the CPF algorithm correctly placed 
the source at a position equidistant from the Hanford and from the Livingston LIGO detectors 90% of the time. It 
was, however, often unable to correctly use the information provided by the Virgo instrument. The CPF algorithm 
can be applied to networks of arbitrary complexity; it is likely that including other interferometers in the network 
will reduce the amount of points on the sky where no instrument can effectively complement the measurements from 
the two LIGO sites. 

Although it has not been discussed in this paper, it should be noted that the CPF algorithm can also be used 
to improve the detection efficiency of a network of interferometers. It is shown in that CPF achieves a higher 
probability of detection at the same false alarm rate than a incoherent approach based only on the comparison of 
event lists from the three interferometers of the HLV network. 
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